m 



PROCEEDINGS 

Of SCIENCE 



MKPH-T-07-18 
SHEP-0744 



O 

o 

(N 

> 

O 

00 



Oh: 



Progress in kaon physics on tiie lattice 



Andreas Juttner*^ 

School of Physics and Astronomy 

University of Southampton 

Southampton, SO 17 IB J 

UK 

E-mail: 



juettnergkph . uni-mainz . de 



> 

en 



CKM-unitarity, direct and indirect CP- violation and the hi ~ \/2 rule in full lattice QCD are 
the focus of this talk. To this end I will discuss and compare recent lattice results for leptonic, 
semi-leptonic and non-leptonic decays of the kaon and neutral kaon mixing and I will motivate 
current best estimates 



o 



X 



//f/A = 1-198(10), /f''(0)= 0.964(5) and B^ = 0.720(39) . 

Moreover new theoretical advances that will improve the quality of these computations in the 
future will be discussed. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 

'New address: Institut fur Kernphysik, Johannes Gutenberg-Universitat Mainz, Johann-Joachim-Becher-Weg 45, 
D-55099 Mainz, Gemany 



© Copyright owned by the ;iuthor(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAHke Licence. 



http://pos.sissa.it/ 



Progress in kaon physics on the lattice Andreas Jiittner 



1. Introduction 

This talk reviews recent results and developments in the phenomenology of the kaon from the 
lattice focussing on simulations of full QCD (i.e. with a degenerate pair of dynamical u and d 
quarks {Nf = 2) and a dynamical 5-quark {Nf = 2 + 1)). Many of the results presented here have 
been computed using different fermion discretisations. This has to be seen as a feature since we are 
entering a phase where in particular observables describing 5J7 (3) -breaking effects in kaons can be 
computed on the lattice with sub-percent level precision. If used to constrain the Standard Model 
and also to constrain models that go beyond, consistency of the results from different discretisations 
will increase the reliability of the lattice output. To this end I want to draw the reader's attention 
to the plenary talks on recent simulations with domain wall fermions (DWF) ^, overlap fermions 
[§], twisted mass fermions (tmQCD) [^] and clover improved Wilson fermions [Q] and also the two 
plenary talks by Creutz ^ and Kronfeld [^ ' discussing different perspectives on the consistency 
of the staggered fermion (KS) formulation. 

Please be aware of the preliminary status of some of the results presented here and that I could 
not cover all recent results and developments of this field. 

This talk is structured as follows: I will start with a discussion of the status of the first-row- 
unitarity of the CKM matrix, particularly focusing on the matrix element \Vus\ which can be deter- 
mined either from leptonic kaon and pion decays or from semi-leptonic kaon decays. I will then 
discuss indirect CP-violation in neutral kaon mixing where the lattice can contribute in terms of the 
computation of the bag parameter Bk- Finally, I wiU discuss direct CP- violation and the A/ = 1/2- 
rule in hadronic kaon decays, where the parameters \£' /£\ and (O can be computed on the lattice. 
While discussing these topics I will also present new theoretical and technical developments which 
will hopefully improve the quality and the precision of these calculations. 

2. The determination of I y„ 5 1 

The three entries Vud, K* and Vub make up the first row of the CKM-matrix [||, ^ and assuming 
its unitarity, |V„dP + |K.vP + iK/jP = 1- Any deviation from 1 on the r.h.s. would be a sign for 
physics beyond the Standard Model. \Vud\ is very well known and neglecting \Vub\ (since it is very 
small) the precision of \Vus\ is currently limiting the accuracy of the unitarity test. 

2.1 I Vui I from leptonic decays 

In 2004 Marciano [ |Io| ] first used the lattice determination of /k/Jk to determine iVusI employ- 
ing the relation 

2 



r{i^Mv„(r)) \v.,P\fJ „„(i-4) 

From the charged kaon and pion life-times and the leptonic partial widths in [ pl| ] one determines 

r{K -^nv^ir)) = 2.528(2) x lO^'^^MeV , (2.2) 

r{K^HV^{Y)) = 3.372(9) x lO^^'^MeV, 



See also Sharpe's plenary talk last year M\. 
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group Nf action 


a/fm 


a from Lm\ nijilMdV 


fKJfn 



QCDSF+ 
UKQCD 

ETM 

MILC 

HPQCD+ 
UKQCD 

NPLQCD 

RBC+ 
UKQCD 

PACS-CS 



[14] 



2 clover (NP) > 0.06 



[[T|, |T6|l 2 max. tmQCD 0.09 
[jnKTim] 2+ 1 Y^^^tc ^ 0.06 

> 0.09 



[20] 



0,1 T^cHISQ 



[21] 



2+1 KSmilc/DWF 0.13 
2+1 DWF 0.11 

2+1 clover (NP) 0.09 



ro 4.2 >300 1.21(3) 

A 3.2 >290 1.227(9)(24) 

A 4 >240 1.197(3)(+t3) 

T 3.8 >250 1.189(7) 

ro 3.7 >290 1.2 18 (2) (+1^) 

Q. 4.6 >330 1.205(18) 

3 >210 1.219(26) 



' lightest pion; NP: non-perturbatively improved; mux. tmQCD: maximally twisted mass QCD 



Table 1: Parameters of gauge configurations from which /k/Zk has been determined. Errors on the results 
are either statistical and systematic or the combined error 



group 



ami -extrapolation 



FVE 



QCDSF+UKQCD 



linear extrapolation 



2 values of a 



ZPT 



ETM 



NLO;t:PT+NNLO (analytic terms) 
and also polynomial 



X^ 



MILC rS;^PT: first NLO + NNLO analytic terms on reduced set 2 volumes 

then include NNNLO-analyt. terms and rS;tPT 

include also larger values of ami, 4 values of a 



HPQCD+UKQCD 



NLO;^PT+NNLO(analytic) 

a^ (conventional+taste breaking)-terms 

3 values of a 



;fPT 



NPLQCD 



NLO;f PT+NNLO analytic terms 



RBC+UKQCD 



NLO;i;PT 



■SU{2)xSU{2) 



2 volumes 



PACS-CS 



NLO;t;PT 



2 volumes, [ 14 ] 



^PT: continuum chiral perturbation theory - rS;fPT: rooted staggered chiral perturbation theory (including cut-off dependence) 

Table 2: Summary of how the systematic effects due to the extrapolation in the quark mass, cut-off effects 
and finite volume effects have been treated in the determination of /k/Zk- 



which then yields 



\^us\ Jk 



:0.07602(23)exp(27)Rc^ 



(2.3) 



with errors from experiment and radiative corrections, respectively. Combining the very precise 
result \Vud\ = 0.97377(ll)exp(15)nuci(19)Rc^ [13] from nuclear j8-decay with the prediction for 
fx/fn from lattice computations one determines |V„.v|. Table |T] summarises the basic parameters 
and results for the gauge field ensembles from which fx/ fn has recently been determined. I begin 
with some comments on the treatment of the systematic error which is summarised in table ^: 

Finite volume effects were studied by all collaborations except NPLQCD by comparing results 
from lattices with different physical volumes but otherwise constant physical parameters (MILC, 
RBC+UKQCD, PACS-CS) and/or using NLO chiral perturbation theory (QCDSF+UKQCD, ETM, 



A preliminary update for this number was given by Marciano at the Kaon 2007 International Conference: |V„j| 
0.97372(10)exp(15)„uci(19)RC Q 
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MILC, HPQCD+UKQCD) [g5|, g. I should comment that the first method is clearly preferred and 
that values of rrij^L w 3 very likely lead to finite volume effects that are not negligible. In addition to 
comparing results from two volumes for some simulation points PACS-CS estimated finite volume 



effects using the approach by Colangelo-Durr-Haefeli []24[]. It is based on Luscher's idea to express 
finite volume effects of pion (kaon) masses in terms of the nn{K) scattering phase shift [ |27| ] and 
predicts larger finite volume effects than NLO ;tPT. Since this approach depends only indirectly 
on chiral perturbation theory through the estimate of the scattering phase shift the predictions are 
expected to be more accurate. 

chiral extrapolation: UKQCD+QCDSF don't see a curvature in their data for fx/ fn and there- 
fore linearly extrapolate to the physical point. Staggered fermions come in tastes and at finite lattice 
spacing the taste symmetry is broken and the corresponding pion spectrum is non-degenerate. The 
continuum chiral effective Lagrangian therefore does not correctly describe the spectrum of stag- 
gered lattice QCD. MILC, who are using AsqTad improved staggered valence and sea quarks, 
therefore employ SU{3) x SU{3) rooted staggered partially quenched chiral perturbation theory 



[ ]28| , [29| ] which simultaneously describes the approach to the chiral and the continuum limit. They 
first fit the NLO expression including NNLO analytic terms to a reduced data set including only 
lighter data points and subsequently extend the fit to the heavier data points by including NNNLO 
analytic terms [p^]. HPQCD-i-UKQCD are simulating the partially quenched theory using HISQ 
[30] valence quarks on the MILC staggered sea [31, 1^]. Arguing that for this setup the taste split- 



ting is reduced sufficiently, they rely on NLO SU(3) x SU(3) continuum chiral perturbation theory 
with added NNLO analytical and cut-off terms. In their recent pubhcation [^] HPQCD-i-UKQCD 
quote the smallest en^or for fx/ fn among all the results presented here and it will be interesting 
to compare the details of their chiral and continuum extrapolation once they are accessible in an 
upcoming publication. ETM, MILC, HPQCD-fUKQCD and NPLQCD also use SU{3) x 5^7 (3) 
NLO continuum chiral perturbation theory adding NNLO analytical terms to the fit-ansatz in order 
to improve the fit quality; PACS-CS is using NLO only. RBC-i-UKQCD use a different strategy. 
Assuming that the strange quark mass is too heavy to be properly described by SU{3) chiral per- 
turbation theory and nix » niu^d they use SU{2) x SU{2) chiral perturbation theory fits for the kaon 
sector [p^], to describe the light quark mass dependence of the decay constants at fixed strange 
quark mass. The results obtained in this way are then interpolated to the physical point of the 
strange quark mass. 

cut-off-effects: Apart from MILC and HPQCD-i-UKQCD who are doing a combined chiral and 
continuum extrapolation only QCDSF-i-UKQCD assess cut-off effects. The latter do not see any 
scaling violations in their data. Some of the remaining collaborations will supplement their current 
results with results on a finer lattice and/or rely on the crude estimates for cut-off effects to be of 
order (aAQCD)^ in the 0(a) -improved theory. 

2.2 \Viis\ from semi-leptonic decays 

The Standard Model expectation for the K -^ nlv (A'/a) semi-leptonic decay rate is [ ^ ] 

^K^nlv = Ci^^/5EW[1 +2A5f;(2) +2Aem]|Ks| V+'^WP . (2.4) 

Here, Cj^ = 1/2(1) is the Clebsch-Gordan coefficient for the neutral (charged) kaon decay and 
/ is the phase space integral which is typically determined from the shape of the experimentally 
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measured form factor p3|]. 5ew is the electro- weak short distance correction and ^su{l) ^^d ^em 
are 5?7(2)-isospin breaking and electromagnetic corrections, respectively. The non-perturbative 
contribution to the process is given in terms of the form factor ff^ (0) defined through the QCD 
matrix element of the vector current V^ = sy^u between the kaon and the pion, 

{n{p„)\V^mK{pK)) = ff''{q^){pK + Pn)^+ f^''{q'){pK- Pn)i,, (2.5) 

where q^ = {pK — p 71)11 is the momentum transfer. Recently the FlaviaNet Kaon Decay Work- 
ing Group determined the very accurate value iKs/f'^CO)! = 0.21673(46) [34| which can only be 



fully appreciated for a determination of 1^^! if lattice computations of f^^{0) reach sub-percent 
precision. 

In chiral perturbation theory ff^ (0) is expanded as 

/f^(0) = l+/2+/4 + ..., (2.6) 

where fx, which corresponds to the leading chiral correction, is fully determined in terms of the 
meson masses [ ^5| ] and takes the value /| ^^ = —0.023 at the physical point. Relevant for lattice 
computations are also the more recent evaluations of /2 in partially quenched chiral perturbation 
theory and the evaluation of finite volume corrections [5^, Evl]. A first estimate of the higher order 



corrections was given in [ ]33| ] and more recently estimates for /4 were given in p8| , p^ , ^(jj \. Thus, 
relying on our knowledge of fx in chiral perturbation theory, on the lattice one computes only the 
corrections, 

Af{mK,mj^) = f+{0,mK,mj^) - {1 + f2{mK,m^)) . (2.7) 

Considering for a moment the first estimate A/ = —0.016(8) for physical pion and kaon masses by 
Leutwyler and Roos in 1984 p3|], it is clear that a precision of 30-40% for A/ is sufficient to reach 
sub-percent accuracy for ff^{0). 

During the last three years a number of collaborations have computed this quantity using 
dynamical fermions. However, in terms of controlling systematic uncertainties the project by the 
RBC-i-UKQCD collaboration which use 2-1-1 flavours of dynamical domain wall fermions with 
pion masses down to 330 MeV is the current state of the art pll,p2, 43] (cf. the summary of recent 



simulations in table g). At this conference also ETM Ill5|, |44|] and QCDSF [^ indicated that they 
will compute the Kit, form factor with twisted mass fermions and clover improved Wilson fermions, 
respectively. QCDSF gave a first estimate in their proceeding [ p5| ] and one has to wait for their final 
error analysis and the extension of their analysis to lighter pion masses and one coarser and one 
finer lattice spacing. 

The technique for the high precision calculation of ff''^ (0) has been set out in [ |3^ and consists 
of a 3-step procedure: 

1) one first computes /^""{qi^^) = /f ''(^Lx) +^max/('«i -'"D/- (^Lx)' where ^^,^ = (m^ - 
ntji)^, from 

{7t\Vo\K){K\Vo\K) {mK + m,f , 2 

{n\Vo\n} {K\Vo\K} ArriKm^ 

The l.h.s. is obtained from ratios of suitable Euclidean three-point functions at large values 
of the Euclidean time. Since both mesons are at rest and due to cancellations of correlations 



and also the cancellation of the vector current renormalisation constants in the ratio in ([2.8|), 
/(f^(^max) '^^i^ t)^ determined with very high precision. 
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group 



Nf 



action a/fm L/fm »%/MeV/^"(0) 



JLQCD 

RBC 

QCDSF 

Fermilab, 
HPQCD, MILC 

RBC+UKQCD 



[46] 
147] 
145] 



clover (NP) 

DWF 
impr. Wilson 



0.09 
0.12 
0.08 



1.8 
2.5 
1.9 



t 



2+1 Wilson li-quark 
impr stag, u Sls 



t t 



[I|,g|,g|]2+1 DWF 0.111.8,2.8 



^550 

>490 

592 

X 
>330 



0.967(6) 

0.968(9)(6) 

0.9647* 

0.962(6)(9) 

0.964(5) 



' computed /q ''(^max) ^"'^ ^h.e.n used slope of form factor from experiment for extrapolation to t/-^ = ("exploratory study"); 

* Error analysis not yet finished 
^ The proceeding in which the result was published does not contain this information 

Table 3: Lattice computations of the Ki^ form factor using dynamical fermions. Errors on the resuhs are 
either statistical and systematic or the combined error. 





^"^-fits 


chiral fits 


0,0 






lin quad pole z 


const lin quad 


-0.2 


^-^W""^"^ 


JLQCD 1 


4^ 


• 


• 


RBC [ 


47] 


o o • 


o • o 


|--0.4 


'""^^^l * 


Fermilab rC^i 
FNAL,MILC ^J 


.+ 


• 


-0.6 




RBC/UKQCDgl] 


o o • o 


o • o 


-0.8 





1.1 . . 1 


i = pole mass froir 


expe 


riment -^ extrapolate to /g (( 


) from/o(9miix) 


0.2 0.4 0.6 0.8 
mj, + ml [GeV] 



Figure 1: Left: Summary of how systematic in chiral extrapolation and ^^-interpolatoin was assessed; right: 
chiral extrapolation of R^f (r.h.s. plot taken from [p3|]). 



2) One then determines f^^iq^) at various other values of q^ < q'^^^ from similar ratios of 3- 
point functions by inducing Fourier momenta (|7?A-,7t| =2n/L, \/^2n/L) into the initial and/or 
final state. In order to interpolate the form factor to ^^ = the collaborations have used first 



or second order polynomial ansatze in q , a pole dominance ansatz /o,poie('? 



/o(0) 



or 



the z-fit (polynomial with improved convergence []49[]) as summarised in the table in fig- 
ure [l|. 
3) In the last step the data for /+'^(0) = fo^{0) has to be extrapolated to the physical value of 
the quark masses and it is common practice to use the ratio 



RAf = 



A/ 



/(f(0)-(l+/2(m^,m;,)) 



{mi- 



ni 



2\2 



mt 



m 



2^2 



(2.9) 



with a polynomial ansatz constant, linear or quadratic in (m^ + m|). Again the table in figure 
[l] summarises which ansatz has been studied by the various collaborations. An example for 
the extrapolation with a linear fit is shown in in the plot in figure [l] which is taken from James 
Zanotti's talk 



Both the ^^-interpolation and the chiral extrapolation are based on phenomenological fit-ansatze 
and are thus sources of systematic uncertainties which can be estimated in terms of the spread of 
results between the different ansatze. The RBC-i-UKQCD-coUaboration [^] has observed a reduc- 
tion in the final error when combining the ^^-interpolation and the chiral extrapolation into one 
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global fit. In the next section we will briefly introduce a new technique 1 5C ] that may entirely re- 
move the uncertainty due to the ^^-interpolation. Before let me make some comments: 

• It would be interesting to have the expression in chiral perturbation theory of the Kij, form 
factor to order /7^ in [39] in terms of the quark masses. This would on the one hand allow to 



understand the slight tension between x?T and the lattice results which I will comment on 
later. On the other hand lattice results could then also be used to constrain or even determine 
the low energy constants appearing in these expressions. 
• Taking the preliminary result by RBC-i-UKQCD [ |4l| ] as the current bench-mark, the contri- 
butions to the overall error on A/ from the lattice decomposes in the following way^ : 

A/ = -0.0161 (45)-^'^'- (15)^ (16)^' (8)'' 
30% 9% 10% 4%, 

for the statistical error, error due to the chiral extrapolation, error due to the ^^-interpolation 
and error due to the finite lattice-cut-off, respectively. In the following we will briefly dis- 
cuss two recent developments which should allow to further reduce the error due to the 
^^-interpolation on the one hand and the statistical error on the other hand. 

2.3 Kb with partially twisted boundary conditions 

Using twisted boundary conditions for the valence quark fields, i.e. q{x + Lk) = e'^''q{x) for 

k = 1,2,3 and leaving the sea quark's boundary conditions untouched (partially twisted boundary 



conditions) [51, 52, p3], one can tune the momenta of hadrons in dynamical simulations in a finite 



box continously. For example the dispersion relation for a charged pion then follows 



^^^ = \jK^ + {p\^t-\0u-ed)Y, (2.11) 

where 0„ and 6d are the twisting angles of the valence up and down quarks, respectively, and p\ai 



is the conventional Fourier momentum. This relation was confirmed numerically in [54]. Twisted 



boundary conditions were also tested in the quenched theory [55] where they have subsequently 
been been applied to study the g^ -dependence of the Kji, form factor using techniques similar to the 



3-step procedure detailed in the previous section [56]. An approach which goes beyond that study 
was developed in [BO] where the matrix element 



{K{pn)\V^mK{PK))=ff''mEK + E^)+f^''{0){EK-E^), (2.12) 

of the 4th component of the vector current is evaluated for two choices of the kinematics, namely 



1) {7imV,\K{dK)%2^,: 10^1 =lJ(4±^)2-,„2, and 4 = 0, 



2) {7i{d,)\V,\K{0)%.^,: \e,\=LJ{'^^f-ml and 0^ = 0. 



T-^ - - (2-13) 



The form factor ff^{0) in eq. (2.12) can directly be extracted from a linear combination of the 
matrix element in the kinematical situations 1) and 2) in (2.13), thus entirely avoiding the q^- 



interpolation and difficult fits to ratios of correlation functions involving the spatial components 



In order to disentangle the various contributions to the systematic error I here quote the results from the individual 
5 -interpolation and chiral extrapolation in Mm. 
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M ~ 540 MeV 

K 

M ~ 300 MeV 



\ ^ 



o point-to-all (fixed source) 
• all-to-all (stochastic source) 



(aq) 



Figure 2: Left: fQ^[q^) with all-to-all propagators; right: the pion form factor fniq^) constructed from 
point-to-all and all-to-all propagators, respectively [H|. 



of the vector current, which are necessary in the 3-step procedure hned out above. The statistical 
errors of the results from the new approach are comparable to the errors using the conventional 
approach. 

2.4 /T/s -decays using propagators from stochastic sources 

Instead of constructing the three point functions relevant to the determination of the ^/3-form 



factor from point-to-all propagators, the ETM collaboration [44] uses all-to-all propagators gen- 
erated with the one-end-trick [|^, 58], thus gaining a volume averaging of the propagator source. 
The pion and kaon momenta in these simulations are induced using partially twisted boundary con- 
ditions [|1|, ^ ^] in combination with the 3-step procedure lined out in section 2.2. Preliminary 
results for this approach were shown in Simula's talk at this conference (cf. l.h.s. of figure ^. Since 
no data for comparison with results from point-source propagators were available for f^^{q^) and 
f^^iq^), one gains some feeling for the improvement by looking at Simula's plot of the pion vector 
form factor f^^{q^) on the r.h.s. of figure ^ 

2.5 Summary for I Ki I 

Figure |3| shows a comparison of recent results for \Vus\ from both leptonic and semi-leptonic 
kaon decays. The grey band in the plot gives the result for \Vus\ which one gets using Leutwyler and 
Roos' result /f''(0) =0.961(8) [B|] and |V„,/f''(0)| = 0.21673(46) by the FlaviaNet Kaon Decay 



Working Group [34]. The lattice results for Ki^ currently support Leutwyler and Roos' prediction. 



Both chiral perturbation theory for /+'^(0) [B9l 59, M and the lattice predictions for /k/ fn tend 



towards smaller values for |V„.v|, thus creating a slight tension. We also see that results for the 
CKM-matrix element from lattice calculations of leptonic kaon decays are becoming competitive 
with results from calculations of semi-leptonic decays. 

RBC-i-UKQCD currently have the best control of systematic effects in the calculation of 
f+^ (0) [|l|, ^ ^, |0|] and I take their result as the current best estimate. 



/f-(0) 



0.964(5) ^|K.| =0.2247(12). 



(2.14) 



Due to the preliminary status of some of the central values and error bars for fx/ fn in table [l| it 
is difficult to compute a world average with a reliably estimated error. As the central value for the 
current best estimate I suggest the weighted average over all results in table |T] assuming Gaussian 
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t^ 



o 

■a 



N^2 _ 




QCDSF+UKQCD 2007 


N^2 - 

N=2+1 
N,=2+1 


■ i 


ETM 2007 

MILC Lattice 2006 


* 




N^2+1 




NPLQCD 2006 


; 


N,=2+1 




HPQCD/UKQCD 2007 


N^2+1 




RBC/UKOCD 2007 (16^) 
RBC/UKOCD 2007 (24^) 

PACS-CS 2007 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


Nj=2+1 






,,,,!,, 


1 1 : 1 1 



N=2 



0.225 



Citigliano et al. 2005 
Jaminetal.2004 
Bijnens & Talavera 2003 

Leutwyler& Roos 1984 

Becirevic el al. 2004 

JLQCD 2005 

RBC 2006 

QCDSF 2007 (no error yet! 

MILC+HPQCD 2004 

RBC/UKOCD 2007 



0.225 



0.23 



0.225 




0.97 



0.975 



Figure 3: Left: Comparison of results for the CKM-matrix element \Vus\ from lattice computations of the 
leptonic and semi-leptonic kaon decays. In the lower plot also results from chiral perturbation theory are 



shown. Right: Compilation of the estimates (2.14) and (2.15) into a combined fit by the FlaviaNet Kaon 
Decay WG (like in [p4|]) to assess whether the first row of the CKM matrix fulfils the unitarity constraint. 



and un-correlated errors^. MILC [|l9|] has carried out the most extensive simulation with a very 
detailed study of the chiral extrapolation, the finite-volume and the cut-off effects all within the 
frame work of rooted staggered chiral perturbation theory and I therefore attach to the averaged 
central value their combined statistical and systematic error, 



h/fn = 1.198(10) -^ |y„,| = 0.2241(24) ., 



(2.15) 



which reveals a slight tension with the experimental value 1.223(12) []11[]. Note that this average 
is dominated by the MILC and HPQCD-i-UKQCD results with staggered fermions which have 
very small errors. The average without the results based on staggered fermions would instead be 
1.211(10). 

The FlaviaNet Kaon Decay Working Group was so kind to provide a version of their unitarity 
fit [34, 61] assuming the above best estimates which is shown in the r.h.s. plot in figure ||^. The 



blue vertical line represents Marciano's update for | V„d| from nuclear j8-decays [12]. The horizontal 
band is the result for \Vus\ from RBC-i-UKQCD's lattice calculation of the semi-leptonic form factor 
[|l|, ^ |6^] and the slightly tilted horizontal band represents (2.15). The solid black line represents 
CKM first-row unitarity (neglecting \Vub\)- This analysis which is partly based on preliminary 
results indicates a tension between the lattice results and CKM-unitarity. With \Vud\ = 0.97372(26) 
and |V„,| = 0.2246(11) one gets \Vud\^ + |K.vP - 1 = 0.0014(7). It will be extremely exciting to 



In the case of assymmetric error bars I have shifted the central value of the result and corrected the error corre- 
spondingly to be symmetric. 

Many thanks to M. Antonelli, F. Mescia and M. Moulson for generating this plot ||62|]. 
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Figure 4: Status of the CKM-triangle by the CKMFitter group Q 

follow the development of this situation as the various collaborations finalize their analysis and in 
particular theoretical errors are further reduced. 

3. CP-violation in kaon systems 

The physical state K^ consists predominantly of the CP-odd K2 and an admixture of the CP- 
even K\ . The decay of the K2 into a CP-even two-pion state is called direct CP- violation and has 
been established experimentally in 1999 by NA48 and KTeV [53, U]. Direct CP- violation which 
occurs when the ^1 decays into the CP-even two-pion state has been established in 1964 by Cronin 



and Fitch []65|]. Direct and indirect CP- violation have been studied on the lattice for many years 
and in the following I will discuss recent advances and results. 



3.1 Indirect CP- violation - neutral kaon mixing 

The experimental measurement of |£/f | = \A{Ki 
a prediction of the kaon bag parameter 



{mz)i=Q)/A{Ks -^ (7171)7=0)1 together with 



B, 



:C(At) 



{K^\Q^=\^i)\K^) 



(3.1) 



define a hyperbola in the plane of the Wolfenstein parameters f) and p which constrains the apex 
of the unitarity triangle shown in figure ^ Here Q'^^^ip.) equals the difference of the parity even 
and parity odd four quark operators Oyy+AA and Oya+^j respectively. Before presenting recent 
results for B^ I will now briefly discuss and compare how the major fermion discretisations fare in 



the computation of the matrix element in eqn. ( |3.1| ). 

Wilson fermions: As a result of the explicit breaking of chiral symmetry due to the Wilson 



term the Bk operator mixes with four other operators O, [pTQ, 



O 



AS=2 



(M) 



'VV+AA 



i^,gl 




VV+AA 



+i:a, 



i}or 



(3.2) 



thus complicating the renormalisation process. It was realised in [|^] that the parity even operator 
Ovv+AA is related to the parity odd operator Oav+va by an axial Ward identity and that the latter 
operator renormalises multiplicatively also for Wilson fermions. While the relevant matrix element 
of the parity even operator can be determined from three-point functions, the relevant matrix el- 
ement for the parity odd operator has to be determined from four point functions. In some sense 
one trades the uncertainties due to the mixing in the case of the parity even operator for noisier 
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signals of the four point functions in the case of the parity odd operator. Both approaches have 



been studied with dynamical fermions in ||69|, [7Q| ]. 



Twisted mass QCD: The the parity even operator in QCD is related to the parity odd operator 
in tmQCD by an axial rotation [^], (^°|Ow+aaI^)qcd = -/(^''|oSVAvl^)tmQCD • Bk can 
now be computed from 3-pt. functions of Ova+av which is expected to give a better signal than the 
4-pt. function which has to be computed in the case of standard Wilson fermions. A benchmark 
computation using this approach in quenched tmQCD has recently been carried out by the ALPHA 
collaboration [^, ^. In their simulation they used degenerate light and strange quark masses, thus 
neglecting possible 56^^(3) -breaking effects. The first case that was studied in that paper was with a 
Wilson 5-quark and two twisted light quarks at twisting angle a = n/2. Although non-degenerate 
s and d quarks are in principle feasible in this approach, here the masses of the s and the d quark 
were tuned to be degenerate, which is non-trivial since they have been discretised differently. The 
second case that ALPHA studied was a = n/A with twisted s and d quarks in which case the 
quarks are automatically degenerate. The simulations were carried out for a number of different 
lattice spacings, thus allowing for a very detailed study of the approach to the continuum limit. The 
authors used the non-perturbatively computed (Schrodinger functional) renormalisation constant 



for OvA+AV [|74|]. One further outcome of this study is that the splitting of the meson spectrum 
due to the explicit breaking of the flavour symmetry with twisted mass fermions reduces as the 
continuum is approached where it is expected to vanish. 



Staggered fermions: Van de Water and Sharpe [75] studied the transformation properties of 
the staggered B/f-operator that couples to external kaons of taste P which correspond to the lattice 
Goldstone kaon. The lattice representation of that operator mixes with many other operators of 



all tastes. In current simulations [J76[ [77|] only the operators with the same taste as the lattice 
Goldstone kaon are actually implemented. The corresponding mixing coefficients are computed 
in perturbation theory. The mixing with other tastes at order a and higher orders in the strong 
coupling constant as well as lattice artefacts entering at order a^ are described by staggered chiral 
perturbation theory: 

^stagg,cont ^ of^^ia?.tQV\ + — [tastcPl +-^ [wrongtastes] + a^Faiitastes] +^2 [all tastes] . (3.3) 

An 4k ^ !^ ^ J- V ^ y 

. ,. unknown 2— loop discretisation 

simulation 

The perturbative coefficients of the terms in a/{47i) are known. Since the higher order perturbative 
coefficients ai^e not known, they are counted conservatively as a^ rather than (a/(47r))^. The above 
expression has 37 free parameters which can be constrained e.g. by measuring the taste-sphtting 
or by first determining a sub-set of parameters at a single lattice spacing. It is also known that 



HYP-smearing [78] reduces the mixing and taste breaking and also improves the convergence 



of perturbation theory and may thus yield a more favourable power counting [79, 8C]. It has to 
be mentioned that non-perturbative renormalisation is in principle possible. However, all current 
simulations of Bk with staggered fermions rely on perturbative renormalization. 

Domain wall fermions and overlap fermions: With chiral fermion formulations the parity- 
even operator Ovv+aa renormalises multiplicatively. For domain wall fermions residual mixing 
with wrong chirality operators was discussed in detail in [pT]]. It is suppressed by (amres)^ and 
therefore negligible (e.g. anises ~ 10^^ for the current RBC-i-UKQCD domain wall fermion data set 
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Nf 



m„L mn/MeV 



_o<.<^- 



NC* 



>°"'ZPT 



Bk 



HPQCD+ 
UKQCD 

Bae, Kim, 
Lee, Sharpe 

RBC+ 
UKQCD 

JLQCD 



[0 2+1 KS^y^c 'iissitf 

2+lKS»lPc 0.125 
2+1 DWF 



0.11 

quenched 
^scalina 

2 gy^ji'y 0.12 




E£|J ^ fix. top 



4.5 


360 


PT 


• 


0.83(18) 


4.5 


360 


FT 


• • NLO 




4.6 


330 


NPR 


• • NLO 


0.720(39) 


2.7 


290 


NPR 


• • NLO 


0.723(12)* 



1" for lightest pioii; ^ statistical error only 



Table 4: Summaiy of current BA:-calculations. 



[|l|]). The chiral symmetry of these actions provides automatic 0(a) -improvement and continuum 
chiral perturbation theory can be used. 

Recent simulations: Current efforts for the calculation of Bk on the lattice with dynami- 
cal fermions are summarised in table ^ HPQCD+UKQCD [76] on the one hand and Bae, Kim, 
Lee and Sharpe pl^ on the other hand are using 2+1 flavour staggered quarks (MILC configura- 



tions [31]) at the same simulation parameters. Both collaborations use HYP smeared [78] valence 



quarks, thus simulating a partially quenched theory. While HPQCD+UKQCD is simulating for 
degenerate s- and J-quarks only, Bae and collaborators also investigate 5'L'(3)-breaking effects 



and fit their data to continuum partially quenched chiral perturbation theory []75|]. It remains to be 
studied in detail whether taste breaking effects for HYP- valence quarks are sufficiently suppressed 
that one can apply continuum chiral perturbation theory instead of staggered chiral perturbation 
theory in order to reliably describe and extrapolate the data. In order to estimate cut-off effects, 
HPQCD+UKQCD compare to simulations at various lattice spacings but otherwise similar simu- 
lation parameters within the quenched approximation. The underlying assumption is that cut-off 
effects in the quenched and unquenched theory behave similarly. RBC+UKQCD [|^, |^, |I|] and 
JLQCD [ p5[ ] use chiral fermions with 2+1 flavours of domain wall quarks and 2 flavours of over- 
lap quarks, respectively. RBC+UKQCD estimate the cut-off effects from the experience with the 
quenched case [81, Q and both collaborations consider non-degenerate s- and J-quarks in the 



partially quenched framework and have renormalised the Bk operator non-perturbatively in the 



RI-MOM scheme []87p. JLQCD's overlap quarks are simulated at fixed topological charge and the 
corresponding finite volume effects are estimated to be at the percent level; this estimate is obtained 
by comparing results from different charge sectors at constant volume and quark mass. JLQCD has 
not yet finalised the error analysis for their value of Bk and it will be interesting to see the the 
full impact of finite volume effects (rrijcL ss 2.7) on the final error budget. Since the results of 
HPQCD+UKQCD and Bae, Kim, Lee and Sharpe have been nicely discussed in last year's plenary 
talk [^ I concentrate here on the ones by RBC+UKQCD and JLQCD. 

RBC+UKQCD simulated on two volumes (L = 1.8fm and L = 2.7fm) and as can be gathered 
from the l.h.s. plot in figure ^ no significant finite volume effects were seen. The plot shows a 
number of partially quenched data points for configurations with sea quark masses that correspond 
to approximately 330 MeV and 420 MeV pions. The data was described and extrapolated using 



NLO SU{2) X SU{2) chiral perturbation theory [32] at fixed values of the strange quark mass which 
was then subsequently interpolated to the physical point (r.h.s. plot in figure ^. JLQCD's results 
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Figure 6: JLQCD results for Bk [BJ 



are illustrated in figure |6[ They used NLO SUi^) x SC/(3) partially quenched chiral perturbation 
theory [75] to extrapolate their data to the physical point. The l.h.s. plot in figure ^ shows a 
significant dependence of the fit-results on the number of data points that were included into the 
fit. After applying the cut mi < nis/l their fits however turned out to be stable and as the rh.s. plot 
of figure ^ suggests that the fits all agree after inclusion of an additional NNLO analytic term into 
their fit-ansatz. 

A summary of all recent lattice computations of Bk with dynamical fermions is given in figure 
^ and I quote the numerical values of the most advanced simulations in table ^. The available results 
for the computation of the B^f -operator Ova+av with Wilson fermions either directly or via the axial 
Ward identity method have rather large errors. The reason is the lack of O (a) -improvement of the 
Bjf-operator, rather heavy light quark masses and perturbative renormalization in these simulations. 
Also the HPQCD-i-UKQCD result from staggered fermions has a large error bar which contains a 
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Bk 

Figure 7: Summary of recent lattice simulations with dynamical fermions; for the JLQCD-result only the 
central value is shown since the error analysis is not yet finished. 



large contribution from the perturbative treatment of the operator mixing. The implementation 
of non-perturbative renormalisation for the staggered B^^-operators should be considered in the 
future. The results with the smallest error bars have all been computed using chiral fermions 
thanks to the multiplicative renormalisation of Oaa+vv which was realized non-perturbatively in 
both simulations. JLQCD currently only quotes their central value for Bk since the error analysis is 
not yet finalized. The current best estimates for Bk from lattice QCD is the one by RBC+UKQCD 



Bk = 0.720(39) , 



(3.4) 



This result is compatible with the result of the previous Nf = 2-simulations with domain wall 
fermions by RBC [^] and also with the central value of the new simulation by JLQCD [^]. 

New developments: Aubin, Laiho and Van de Water [ p^ ] developed the partially quenched 
chiral perturbation theory for domain wall valence fermions combined with AsqTad staggered sea 
quarks. In this mixed action ansatz [ pO| ] the symmetry properties of the domain wall valence quarks 
protects the B^^-operator from mixing with operators of non-trivial taste structure. Compared to 



continuum partially quenched chiral perturbation theory [ ]75| ] there are only two additional param 



eters. With the large set of staggered sea quark configurations by the MILC collaboration [ |19| ] 
a computation of Bk using this approach seems feasible and very preliminary results have been 



presented at this conference []91p. 

The scale evolution of the Bk operator is usually carried out in perturbation theory. In order 
to remove the uncertainty due to perturbation theory the ALPHA collaboration has determined 
the scale evolution of the parity-odd operator Ova+av non-perturbatively in the continuum limit 
of the Nf = 0,2 Schrodinger functional scheme [32, ^]. If Ova+av is also renormalised in the 
Schrodinger functional scheme, the discretisation independent result for the scale evolution can 
be used to compute the running of Bk as obtained with Wilson fermions using the axial Ward 
identity method and the running of Bk as obtained using twisted mass QCD or Ginsparg-Wilson- 
type fermions. 

Assuming that Ova+av in the continuum limit of twisted mass QCD has been non-perturbatively 
renormalised at the scale /i, it was suggested in [ p4| ] that the renormalisation constant for the corre- 
sponding operator determined using Ginsparg- Wilson-type fermions, Z^^ {pL^go), at non- vanishing 



14 



Progress in kaon physics on the lattice Andreas Jiittner 

lattice spacing could be defined via 

OvA+Av(M,mp5)|c.l.oftmQCD = Z*^^(At,^o)Oy^^y(go,'"P5) + 0{a^) , (3.5) 

once mps on the l.h.s. and r.h.s. have been tuned to the same value. The identity holds due to the 
automatic 0{a) improvement of chiral fermions. 

3.2 The Holy Grail - direct CP-violation 

CP- violation is described in terms of the isospin amplitudes 

AiK'' ^ n+n-) = sj^A^e'^ + ^1^2^'^ and A(/:0 ^ n^n"") = sj^A^e'^ - ^^Aje'* (3.6) 

from which one constructs the parameters e' = -J^e"^ ( ^^ — ^^ j and (0 = ||^ describing direct 
CP-violation and the A/ = 1 /2-rule, respectively. The isospin amplitudes are defined in terms of 
the matrix element 

(7171(7)1 - iJ^\K^) = A/e'* . (3.7) 

where the relevant effective Hamiltonians ||9^, Q are 

%KrfV;vLC,-(M)0/(M), 
^^ (=1 



f7=± 

10 (3.8) 



for the four- and three flavour case, respectively, where G^j ^^^^ ^1,2 ^^ current-current opera- 
tors, 03,4,5,6 are QCD penguin operators and 07,8,9,10 are EW penguin operators. .^^^^ does not 



contain penguin diagrams. Giusti et al. [ ]970 are studying this Hamiltonian in order to qualitatively 



assess the role of the charm quark in the A/ = 1/2 rule [97] and their programme has been nicely 
reported in Hernandez's plenary talk at Lattice 2006 []98|]. Large scale dynamical lattice simula- 
tions with chiral fermions with the aim to compute the physical values of le'/^l ^^i^ ^ ^^^ ^^^"' 
where the physical charm quark that is too heavy for current lattice simulations has been integrated 



out. Two-pion final states are notoriously difficult to handle on the lattice [99]. Instead of com- 
puting (7171(7)1 — iM'\K^) directly one uses chiral perturbation theory at LO and NLO in order to 
relate the matrix elements of interest Xa K ^ vacuum, K ^ % and K ^ %% matrix elements at 



unphysical kinematics which can be computed on the lattice more easily [|100|, 101,|102, 103]. In 



two impressive works by the CP-PACS [ |104| ] and the RBC collaboration [ ]105| ] the approach has 
been demonstrated to work. The use of the quenched approximation and the rather large values 
of the light quark masses in these calculations of course have to be overcome. To this end Bae, 

(3/2) 

Kim and Lee carried out exploratory studies of the calculation of O7 g using A/y = 2 + 1 flavours 



of staggered quarks [ 106[ , |107| ]. Very recently the RBC-coUaboration has started to repeat their 



calculation but this time with A'/ = 2 + 1 flavours of dynamical domain wall fermions. Their lattice 
spacing is a~' = 1.73GeV and the spatial size of the lattice is L = 2.7fm. The collaboration is 
planning to simulate for 2 sea quark masses corresponding to pion masses of m-ji w 330MeV and 
410MeV, combined with 6 partially quenched valence quark simulation points and a fixed strange 



quark mass. All operators will be renormalised non-perturbatively (RI/MOM scheme [87]). One 
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1 matrix elements: First row: Signal for 
and TniTT/Yrpy ^nd the mass dependence of the latter. Second row: The power subtrac- 



of the many technical difficulties that arise in the project is the appearance of power divergencies 
which have to be subtracted like e.g. 

{n+\0f"-^\K+) + r]u{ms-m,m{sd)\K+), 



(3.9) 



where the subtraction coefficient 77 1, is determined from the slope of the quark mass dependence 
of the ratio of matrix elements 

{0\Oi\K^) 



VO,i + T]l„ 



m, 



■md). 



(3.10) 



The results for the involved matrix elements and the mass dependence of the ratio are illustrated 
in the plots in the first line of figure ||. That the subtraction ( |3. 10| ) works very nicely in practice 
is illustrated in the second line of figure p. The residual chiral symmetry breaking of domain 
wall fermions leads to an ambiguity in the normalisation of the subtraction procedure. This is 
however irrelevant here since only the slope of the subtracted matrix element with respect to the 
quark masses is of interest. The example given here is representative for the technical challenges 
encountered in this project which the RBC-coUaboration is now carrying over to the NLO-level of 
chiral perturbation theory. 

For completeness I also want to mention the study in [108] where the K ^ nn matrix element 



,(3/2) 



for Oj n was evaluated directly for unphysical kinematics and also another exploratory study of 



K ^ TiK at NLO using Nf = 2 flavours of domain wall fermions in [ |109| ]. 

Another technical difficulty in the calculation of K ^ nn processes in the A/ = 1/2-channel 
is the computation of quark-disconnected diagrams that appear as a result of Wick contractions. 
Kim and Sachrajda [110] propose to compute the process Kn^ -^ n^ instead of K ^ n^n^. It 
turns out that the time-like disconnectedness in the latter process is replaced by a space-like dis- 
connectedness in the former which is technically much easier to handle. The relation between the 
two processes has been worked out in NLO chiral perturbation theory. 
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4. Summary 

In this talk I have reviewed recent progress in kaon physics. I discussed the determination 
of |V„.«| from lattice predictions for JkI fn and for ff^{0). The computation of the latter with a 
precision of 0.5% in full QCD is now feasible. The use of all-to-all propagators in the calculation 
of the relevant meson correlation functions and a new approach using partially twisted boundary 
conditions will very likely further reduce the error in the near future - new calculations based on 
these techniques are under way. The determination of \Vus\ from fx/ fn will become an indepen- 
dent competitor for the determination from ff^{Q). I expect the current precision of the various 
results (the best being 0.6% in the case of HPQCD/UKQCD) to be updated in the near future as 
collaborations continue extending their sets of gauge field configurations towards lighter masses 
and also towards smaller lattice spacings. This will then also allow to further improve the quality 
of the combined analysis of \Vus\ from /k/ fn and from ff^{0) in view of CKM-unitarity which 
will be exciting to monitor. 

For Bk the results from the various collaborations agree within errors. The magnitude of 
the errors however varies strongly and the results from chiral fermion formulations currently look 
most promising with an error of about 7% in full QCD. Collaborations using staggered quarks in 
their calculation of Bk should implement non-perturbative renormalisation in order to reduce the 
currently rather large error bars. 

For AS' = 1 the RBC collaboration is in the middle of a large scale project with light domain 
wall fermions in large volume and it will be interesting to see whether they can reproduce the mea- 
sured value of Is'/^l and shed light on the A/ = 1/2-rule - the simulations are however technically 
extremely demanding. 
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